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Gross-Pitaevskii equation with gain is used to model Bose Einstein condensation (BEG) fed by the 
surrounding thermal cloud. It is shown that the number of atoms continuously injected into BEG 
, from the reservoir can be controlled by applying the external magnetic field via Feshbach resonance. 
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^ ■ I. INTRODUCTION 

' Gross-Pitaevskii (GP) equation is one of the most important nonlinear equations used in physics. It appears in 
many branches of physics such as nonhnear optics, fluid mechanics, condensed matter and theory of superfluidity. 
I— —I For example, GP equation is used to explore the macroscopic behavior of the Bose-Einstein condensation (BEG) at 
^ low temperature. Its validity in the BEG is based on the condition that the s-wave scattering length should be much 
smaller than the average distance between the atoms and that the number of atoms should be much larger than 
I one. Experimental achievements and theoretical studies of Bose-Einstein condensates (BEG) of weakly interacting 
. ' atoms have stimulated intensive interest in the field of atomic matter waves and nonlinear excitations such as dark 
"c^ ! QIMI^ and bright solitons 4, dIBiH- The dynamics of BEG crucially depends on the sign of interatomic interaction, 
i: Dark (bright) solitons that are solutions to the GP equation can be generated with repulsive (attractive) interactions, 
I ' resulting from the positive (negative) s-wave scattering length. In the present study, we start by considering the 
. Gross-Pitaevskii (GP) equation modified by a linear term. The Gross-Pitaevskii gain equation reads; 

a ■ 

8 ■ i^^ = -|^V^* + V,ra,^ + ^a,|v|/p* + ^n^^ , (1) 
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^-H \ where Vtrap — — uo^lx^ + y^) H ^1^"^ i o.s{t) is the time-dependent scattering length and the constant F is the gain 

^ \ (loss) term for F > (F < 0). The complex linear term with positive F accounts for the mechanism of loading 
r the thermal cloud into the BEG by optically pumping atoms from the external cold atomic-beam source. From the 
\^ physical point of view, replenishing is required in steady state to compensate for various intrinsic loss mechanisms, 
. such as collisions with hot atoms from the background vapor and inelastic two and three-body collisions between 
' the trapped atoms. However, for the negative values of F, the above equation describes a BEG that is continuously 
, depleted by loss. 

' For attractive interaction, the number of atoms in the condensate grows until the total number of atoms in the 
condensate exceeds a critical value with which the BEG undergoes collapse. During the collapse, the number of 
atoms is decreased. If the condensate is fed by a surrounding thermal cloud, then the condensate undergoes cycles 
of growth and collapses. Hulet and his team at Rice university observed the growth of a condensate of trapped ^Fi 
. atoms with attractive interaction and its subsequent collapse 8] . Ketterle with his team measured the rate of growth 
^ of a ^•^Na condensate fed by a thermal clo ud 131 . Growth of a Bose-Einstein condensate from thermal vapor was also 



Q , experimentally realized for atoms of *^Rb [10|. Some methods have been introduced theoretically to account for the 
O ■ growth of the BEG. GP gain equation was used to model the growth of a BEG by Drummond and Kheruntsyan Tlj. 
^ ' They determined that, as the condensate grows, the center of mass oscillates in the trap. 

Furthermore, GP gain equation is of special importance in the field of an atom laser, that is a device which produces 
an intense coherent beam of atoms by a stimulated process %2, 13]. 

In recent years, GP gain equation has also been studied mathematically by some authors [13. Iislll^ . GP gain 
equation can also be used in the theory of optical fibers if the harmonic trap potential is assumed to be zero [islll^. 
The interaction effect in the system is determined by the s-wave scattering length. In ref . ^3 ; it has been shown 
that controlling the generation of bright and dark soliton trains from periodic waves can be achieved by the variation 
of the scattering length. Recent experiments have demonstrated that tuning of the s-wave scattering length can be 
achieved due to the Feshbach resonance It offers a possibility to vary the interaction strength in 

ultracold atomic gases simply by applying an external magnetic field. 

as{t)^a{l + ^/[B^-B{t)]} , (2) 

where a is the asymptotic value of the scattering length far from resonance. Bit) is the time-dependent externally 
applied magnetic field, A is the width of resonance and Bq is the resonant value of the magnetic field. The strength 
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of the interaction can be adjusted experimentally from large negative values to large positive values. As a result, the 
experiments with magnetically trapped ultra cold atomic gases, where the s-wave scattering length fully determines 
the interaction effects, have an unprecedented high level of control over the interatomic interactions. With this 
experimental degree of freedom, it is possible to intensively study the dynamics of BEC. 

Here, we will study the changes in the number of atoms supplied by a surrounding thermal cloud due to the variations 
of the scattering length by the external magnetic field. We will assume that F does not change but the s-wave 
scattering length changes via Feshbach resonance. 

By supposing the number of particles is not very large for the cigar-shaped traps, it is legitimate to use the reduced 
ID GP equation 
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where = lo\Ilo\, F — ^ /uo±_ and is the Bohr radius. Here, the coordinate z and the time t are measured in units 
of a±_ and \/lo^, respectively. The relationship between the macroscopic wave functions ^'(r, t) and "^{z,t) is given 

by 



*(f,t) = ^=^^exp(-ic^ii-^-^)x*(— , (4) 



where a± = {h,/jnuj±Y^'^ . To extract qualitative information and to improve the understanding of the underlying 
physics, we will transform the GP gain equation Q analytically to the one with the constant scattering length. 
To get rid of the time-dependent scattering length from the equation the two transformations arc introduced as 
follows; 

^ f{t)-^'^ cxp 1^ j ^^^^7 ^'^ dt' + ^ z^^$(z,t); z~.f{t)z, (5) 
where dot denotes time derivation and the time-dependent function f{t) is given by 
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where $7 is a constant. Note that the time-dependent function f{t) is chosen in such a way that it satisfies 
f^f + Lu^f'^ = 0^. The scale transformation on the coordinate accounts for the contraction or the expansion of 
the condensate depending on the behavior of f{t). The time-dependent scaling of the coordinate and re-scaling of 
the wave function have been previously applied to harmonic oscillator potential with the time-dependent angular 
frequency (t) jl^, 0, ll^ . These scale transformations are used here to explore the effect of Feshbach resonance on 
the number of atoms injected into BEC. 

Now, let us first substitute the transformation for \1/ and then for the coordinate into the equation Note that 
under the coordinate scaling, the time derivative operator transforms as d/dt^ d/dt— {f/f)z d/dz. We observe 
that the equation Q can be reduced to the another GP gain equation with the constant scattering length 
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where a is a constant scattering length. We made another transformation on time t — f dt in 

. -^9 . 

the last step. The necessary and sufficient condition for the existence of such a reduction is given by 
as{t) = a exp ^—2 J "^^'^^^ — ^dt'^. As a special case, if >> 1 for some special values of ut'^ and il^, then 
as{t) can be approximated as as{t) « a f~^e~'^^*. 

The magnitudes of trapping frequencies [fl, lu) are scaled. We have freedom to make them equal by choosing Q ~ uj 
or to turn-off the trap by choosing $1 = 0. 

It is interesting to observe that the equation (PJ can be transformed into the another GP gain equation with a differ- 
ent time parameter ((JJ. The two equations (|3I7(I describe the two different condensates continuously replenished by 
pumping from the reservoir. The only physical difference between the two condensates is that the latter one is with 
the constant scattering length while the former one is with Feshbach resonance. The transition from one to the other 
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is given by the transformations lO- 

This reduction is of great importance in the sense that ahhough the corresponding exact solutions are not known, it 
allows us to understand the effect of Feshbach resonance on the dynamics of the condensates. 

Let us find how the Feshbach resonance changes the number of atoms supplied to the system from the reservoir by 
using jsj. Note that the gain coefficient 7 doesn't change for the two condensates. 

1*12^3^ ^ exp f / 2Ml_2ldA [ m^d\ . (8) 



p 

This formula specifically shows the effect of the time-dependent nature of the scattering length on the number of 
atoms in the condensates. The relation © has been derived without obtaining the exact solution to the GP gain 
equation. It is concluded that simply by applying the external magnetic field, one can control the number of atoms 
transferred from the surrounding cloud. In an experiment for the negative sign of 7, one can also take advantage of 
Feshbach resonance to control the number of atoms left in the BEG. 

The method applied here is analytical and only works for a special as(t). Numerical calculations can be performed 
for the most general as{t) on the basis of the fact that Feshbach resonance can be used to control the number of 
atoms supplied to the system. 

As a special case, if we set 7 = 0, then the number of atoms doesn't change for the two BEGs 13I7|I as expected. In 
other words, Feshbach resonance does not change the number of atoms if there is no reservoir. The transformations 
(O conserve the total number of atoms. However, the presence of f{t) in the transformations Q has the effect of 
increasing the central peak of the condensate with a consequent contraction of the atomic cloud. 
If 7 7^ 0, the number of particles as well as the condensate size is changed by Feshbach resonance as can be seen from 
the equation ijHJi. The transformations (5) allow us to understand the underlying physical effects of time-dependent 
nature of the scattering length on GP gain equation, which is used to model the condensate growth and atom laser. 
Let us study the evolution of the system if we choose = 0. In this case, the equation {TJ describes a BEG with 
the trap potential turned-off and continuously replenished by pumping from the reservoir. However, the equation © 
describes another BEG whose scattering length changes with time through the application of the external magnetic 
field. Atoms are injected into both condensates by the same gain coefficient 7. Since the harmonic potential does 
not exist and the scattering length is constant in Q, its solution can be constructed in terms of the well-known 
bright and dark solitons depending on the sign of the scattering length. Then, the exact solution for Q can also be 
constructed by using the analytical transformations given in (Q. 

The time parameter transforms as r = -; ^ when il = 0. The transformations (jSJ works in the time 

sincjt -|- cos Lut 



interval < ujt < 37r/4. In this interval, f{t) — ^(f -I- sin2a;i)/w varies from y/l/uo to zero. Switching-off the trap 
potential has dramatic effects on the dynamics of the system. The longitudinal width of the two condensates can be 
compared by making use of the coordinate transformation (z — > z \/ {1 ~\- sin 2Lot) jio) . If it is equal to L for the BEG 
described by Q, then it is equal to L y^T+sin^wi)/^ for the other BEG described by (PJ. The time = 7r/4 
is a turning point in the sense that the longitudinal width of the BEG starts to decrease from this time while it is 
increasing in the interval Q <u}t < -k /A. 

The condensate collapses into itself at time uit — 'in/A since the longitudinal width of the BEG © goes to zero. 
Furthermore, the density of the BEG shows asymptotic behavior depending on the sign of 7. At time ojt = 37r/4, 
l^'P oo|$p for 7 > while the longitudinal width of the BEG Q approaches zero. The matter wave function 
becomes somewhat like a Dirac delta function. This is reasonable, since atoms are constantly injected into the 
condensate. However, this is not physical, because the BEG becomes a singular point. From the physical point 
of view, GP gain equation is not valid for the description of the BEG if the density of the BEG is big enough. 
To prevent the formation from being a true singularity, the three-body recombination losses are expected to cause 
rapid depletion of the BEG. In order to model atom loss due to 3-body recombination, a phenomenological term 
proportional to 4"^ should be added to the GP gain equation (1) for a good physical description of the BEG. Note 
also that, the scattering length as{t) goes to zero if 7 > at this time. Physically, the scattering length can be made 
zero by tuning the external magnetic field. For example, at a magnetic field oi B k, 166.5 G, the scattering length 
vanishes for ^^Rb and the gas behaves effectively as an ideal Bose gas. However, for the negative values (7 < 0), 
|^|2 _^ o|$|2 at time ujt = 37r/4 . Since the atoms are lost for the negative values of 7, both the longitudinal width 
and the number of atoms get smaller in time. The scattering length as{t) goes to infinity if 7 < at this time. 
Let us now study the case if > 0. The equations I|3I7|I describe two condensates in the confining harmonic traps 
with two different trap frequencies a;^, VP' . The scattering length is time-dependent for the former, while it is 
constant for the latter. Both of them arc continuously replenished by pumping from the reservoir with the same gain 
coefficient. 

The relation between the time parameters is given by VLlanVLT = f + ^/f + fi^ tanojt, where Q < tot < 7r/2. /(t) 
is increasing in the interval < i < 7r/2 and decreasing in the interval 7r/2 < t < tt/4. In this case, f{t) does not 
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become zero. Then, the singularity in the BEC doesn't occurs in this case. This is because of the existence of the 
harmonic trapping potentiaL The harmonic potential which has zero point energy in the ground state plays a central 
role in preventing the condensate from being singularity. 

In conclusion, although the atoms are loaded with the same gain coefficients, the time dependent nature of the 
scattering length changes the number of atoms transferred to the BEC from the reservoir as can be seen from (8). 
The transformations introduced here provides for an analytic way to investigate the effect of Feshbach resonance on 
the number of atoms which are continuously loaded or depleted by loss. 
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